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(Dated: June 3, 2015) 

Using data samples collected with the BESIII detector at the BEPCII collider, we measure the 
Born cross section of e + e” —> pp at 12 center-of-mass energies from 2232.4 to 3671.0 MeV. The 
corresponding effective electromagnetic form factor of the proton is deduced under the assumption 
that the electric and magnetic form factors are equal (|Ge| = \Gm\). In addition, the ratio of electric 
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to magnetic form factors, \Ge/Gm\, and \Gm\ are extracted by fitting the polar angle distribution 
of the proton for the data samples with larger statistics, namely at y/s = 2232.4 and 2400.0 MeV 
and a combined sample at yfs = 3050.0, 3060.0 and 3080.0 MeV, respectively. The measured cross 
sections are in agreement with recent results from BaBar, improving the overall uncertainty by 
about 30%. The \Ge/Gm\ ratios are close to unity and consistent with BaBar results in the same 
q 2 region, which indicates the data are consistent with the assumption that |Gb| = \Gm\ within 
uncertainties. 

PACS numbers: 13.66.Be, 14.20.Dh, 13.40.Gp 


I. INTRODUCTION 

Electromagnetic form factors (FFs) of the nucleon pro¬ 
vide fundamental information about its internal struc¬ 
ture and dynamics. They constitute a rigorous test of 
non-perturbative QCD as well as of phenomenological 
models. 

Proton FFs can be measured in different kinematic re¬ 
gions by i) lepton-proton elastic scattering (space-likc, la¬ 
beled SL) ii) electron-positron annihilation into a proton- 
antiproton pair or proton-antiproton annihilation into 
an electron-positron (time-like, labeled TL). The low¬ 
est order Feynman diagram of lepton-proton scattering 
is shown in Fig. 1(a). The momentum transfer squared, 
q 2 , is negative and the FFs are real functions of q 2 . 
The lowest order e + e _ annihilation process is shown in 
Fig. 1(b). Here, q 2 is positive and the FFs are complex 
functions of q 2 . The basic kinematic variables are also 
shown in Fig. 1, where fc, k' are the initial and final elec¬ 
tron momenta and p , p' are the initial and final proton 
momenta. Since the electromagnetic vertex of the lepton 
is well known, one can reliably extract the proton electro¬ 
magnetic vertex by measuring the cross section and 
the polarization. Assuming one-photon exchange, i.e. in 
Born approximation, and under the basic requirements 
of Lorentz invariance, the hadronic vertex can be param¬ 
eterized in terms of two FFs, F± and F2 [1], 

r m (p',p) = 7/^i O? 2 ) + Me 2 ). W 

where the element < 7 pl , = 7 M 7„ — 7^7^ is a representa¬ 
tion of the Lorentz group, m p is the mass of the proton, 
k p = 9p 2 2 is the anomalous magnetic moment of the 
proton, g p = fy-, p p = 2.79 is the magnetic moment of 
the proton and J = 1 is its spin. The functions F\ and 
F2 are the so called Dirac and Pauli FFs, respectively. 
The optical theorem, applied to lepton-nucleon scatter¬ 
ing, implies that at the lowest order the FFs are real in 
the SL region [2] [3], i.e. the complex conjugate of the 
amplitude in Fig. 1(a), A4^, is identical to Ai. In the TL 
region, as in in Fig. 1(b), the FFs can be complex above 
the first hadronic threshold, that is, above twice the pion 
mass. 

The Sachs FFs, electric Ge and magnetic Gm, are 
introduced as linear combinations of the Dirac and Pauli 
FFs [4]. Concerning the SL region in the Breit frame, Ge 


and Gm are the Fourier transforms of the charge and 
magnetization distribution of the nucleon, respectively. 
Gm and Ge are proportional to spin-flip and non spin- 
flip amplitudes, respectively. They are expressed as 

G E (q 2 ) = F 1 (q 2 ) + ^ Kp F2(q 2 ), ( 2 ) 

Gm(? 2 ) = Fj(g 2 ) + K p F 2 (q 2 ). (3) 

In the TL region, the center-of-mass (c.m.) system is 
equivalent to the Breit frame since the helieities of 
baryons are opposite for the spinors aligned in Gm and 
are the same for the spinors aligned in Ge- 

In the SL region, FFs have been extracted by the 
Rosenbluth separation method [5], as well as, more 
recently, by the recoil proton polarization transfer 
method [ 6 ]. The latter has been applied to obtain the 
PpGe/Gm ratio. Results from the GEp-II experiment 
at JLab’s Hall A [7, 8 ] for p p Ge/Gm show that this ra¬ 
tio decreases rather quickly with increasing Q 2 , where 
Q 2 = —q 2 > 0, while results achieved by the Rosenbluth 
method show an almost constant ratio [9]. The discrep¬ 
ancy between the Rosenbluth and the polarization trans¬ 
fer method may be resolved by including higher order 
corrections like two-photon exchange. A small correction 
to the Rosenbluth separation could imply a large correc¬ 
tion for the extraction of Ge, since Ge is the slope of 
the Rosenbluth plot. However, the correction of includ¬ 
ing two-photon exchange is small and cannot significantly 
influence the results of the polarization transfer experi¬ 
ment. 

In the TL region, measurements have been performed 
in the direct production channel e + e _ —> pp [10-14], in 
the radiative return channel e + e _ —► pp^isr) [15, 16] 
where 7 isr refers to a photon emitted by initial state 
radiation (ISR), and in pp —> e + e~ [17-19] experiments. 
In cases where the data sample is too small to extract an¬ 
gular distributions and disentangle \Ge | and |Gm|, the 
effective proton FF |G| can be calculated from the to¬ 
tal cross section, assuming \Ge\ = \Gm\- This assump¬ 
tion is valid at the pp mass threshold, if analyticity of 
the FFs holds, implying that at threshold the angular 
distribution should be isotropic. In the PS170 exper¬ 
iment at LEAR [17], the effective proton FF was ob¬ 
tained, as well as the \Ge/Gm\ ratio, from pp threshold 
up to y/s = 2.05 GeV. In the BaBar experiment at PEP- 
II [15, 16], the cross section was measured using the ISR 
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(a) (b) 




FIG. 1. (a) Feynman diagram of ep —» ep elastic scattering at the lowest order, (b) Feynman diagram of 
e + e“ —> pp annihilation at the lowest order (identical to that of the reverse reaction pp —> e + e~ with e <A p 
exchange.) 


method from the pp production threshold up to y/s = 6.5 
GeV. The \Ge/Gm\ ratio was measured from threshold 
up to y/s = 3.0 GeV and the result shows an inconsis¬ 
tency with respect to the PS170 results. 


The presence of vector resonances, like p, u> and </> in 
the unphysical region , below the pp threshold, can in¬ 
fluence the functional form of the FFs in the physical 
region. Hence the FFs, in particular the ratio \Ge/Gm\, 
in the TL region cannot be simply extrapolated from the 
SL ones. Until now it has been assumed that all FFs 
respect analyticity, which should allow to calculate their 
behavior in the unphysical region thanks to dispersion re¬ 
lations [20] using the available data in both the TL and 
SL regions. In the SL region, the ratio p p Ge/Gm has 
been measured at 16 Q 2 values in (0.5, 8.5) GeV 2 with 
the best precision to 1.7% [7, 8], while the present pre¬ 
cision of |G_e/Gm| in the TL region exceeds 10% by far. 
Therefore, it is necessary to improve the measurement of 
\G e /G m \ ratio in the TL region. 


The experimental determinations of proton FFs are 
important input for various QCD-based theoretical mod¬ 
els. There are plenty of theoretical approaches applied to 
explain TL FFs: Chiral Perturbation Theory [21], Lattice 
QCD [22] [23], Vector Meson Dominance (VMD) [24], the 
Relativistic Constituent Quark Model (CQM) [25], and, 
at high energies, perturbative QCD predictions [26]. 


In this paper, we present an investigation of the process 
e + e _ —> pp based on data samples collected with the Bei¬ 
jing Spectrometer III (BESIII) [27] at the Beijing Elec¬ 
tron Positron Collider II (BEPCII) at 12 c.m. energies 
{y/s). The Born cross section at these energy points are 
measured and the corresponding effective FFs are deter¬ 
mined. The ratio of electric to magnetic FFs, \Ge/Gm\, 
and \G m \ are measured at those c.m. energies where the 
statistics are large enough. The results are consistent 
with those from BaBar in the same q 2 region. 


II. THE BESIII EXPERIMENT AND DATA 
SETS 

BEPCII is a double-ring e + e _ collider running at 
c.m. energies between 2.0-4.6 GeV and reached a peak 
luminosity of 0.85 x 10 33 cm~ 2 s _1 at a c.m. energy of 
3770 MeV. The cylindrical BESIII detector has an effec¬ 
tive geometrical acceptance of 93% of 47 t and is divided 
into a barrel section and two endcaps. It contains a small 
cell, helium-based (40% He, 60% CaHg) main drift cham¬ 
ber (MDC) which provides momentum measurement for 
charged particles with a resolution of 0.5% at a momen¬ 
tum of 1 GeV/c in a magnetic held of 1 Tesla. The 
energy loss measurement ( dE/dx ) provided by the MDC 
has a resolution better than 6%. A time-of-flight system 
(TOF) consisting of 5-cm-thick plastic scintillators can 
measure the flight time of charged particles with a time 
resolution of 80 ps in the barrel and 110 ps in the end- 
caps. An electromagnetic calorimeter (EMC) consisting 
of 6240 Csl (Tl) in a cylindrical structure and two end- 
caps is used to measure the energies of photons and elec¬ 
trons. The energy resolution of the EMC is 2.5% in the 
barrel and 5.0% in the end-caps for photon/electron of 1 
GeV energy. The position resolution of the EMC is 6 mm 
in the barrel and 9 mm in the end caps. A muon system 
(MUC) consisting of about 1000 m 2 of Resistive Plate 
Chambers (RPC) is used to identify muons and provides 
a spatial resolution better than 2 cm. 

Monte Carlo (MC) simulated signal and background 
samples are used to optimize the event selection crite¬ 
ria, estimate the background contamination and evalu¬ 
ate the selection efficiencies. The MC samples are gen¬ 
erated using a GEANT4-based [28] simulation software 
package BESIII Object Oriented Simulation Tool 
(BOOST) [29], which includes the description of geom¬ 
etry and material, the detector response and the dig¬ 
itization model, as well as a database of the detector 
running conditions and performances. In this analysis, 
the generator software package Conexc [30] is used to 
simulate the signal MC samples e + e~ -A pp, and cal¬ 
culate the corresponding correction factors for higher 
order process with one radiative photon in the final 
states. Another generator Phokhara [31] serves as a 
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cross check of the radiative correction factors. At each 
c.m. energy, a large signal MC sample with more than 
10 times of the produced events in data for the process 
e + e _ —> pp, contributing 0.15% statistical uncertainty 
on the detection efficiency, is generated. Simulated sam¬ 
ples of the QED background processes e + e _ -A l + l~ (1 
= e, p) and e + e _ — > 77 are generated with the genera¬ 
tor Babayaga [32], The other background MC samples 
for the processes with the hadronic final states e + e _ —> 
h + h~ (h = 7 r, K), e + e~ -A ppir 0 , e + e~ -A ppTT°n° and 
e + e _ —> A A are generated with uniform phase space dis¬ 
tributions. The background samples are generated with 
equivalent luminosities at least as large as the data sam¬ 
ples. 


III. ANALYSIS STRATEGY 

A. Event selection 

Charged tracks are reconstructed with the hit informa¬ 
tion from the MDC. A good charged track must be within 
the MDC coverage, | cos 6\ < 0.93, and is required to pass 
within 1 cm of the e + e~ interaction point (IP) in the 
plane perpendicular to the beam and within ±10 cm in 
the direction along the beam. The combined information 
of dE/dx and TOF is used to calculate particle identifi¬ 
cation (PID) probabilities for the pion, kaon and proton 
hypothesis, respectively, and the particle type with the 
highest probability is assigned to the track. In this anal¬ 
ysis, exactly two good charged tracks, one proton and 
one antiproton, are required. To suppress Bhabha back¬ 
ground events, the ratio E/p of each proton candidate is 
required to be smaller than 0.5, where E and p are the en¬ 
ergy deposited in the EMC and the momentum measured 
in the MDC, respectively. The cosmic ray background is 
rejected by requiring |T trk i - T trk2 | <4 ns, where T trk i 
and Ttrk 2 are the measured time of flight in the TOF 
detector for the two tracks. For the samples with c.m. 
energy /s > 2400.0 MeV, the proton is further required 
to satisfy cos 6 < 0.8 to suppress Bhabha background. 

After applying the above selection criteria, the distri¬ 
butions of the opening angle between proton and an¬ 
tiproton, Opp , at c.m. energies /s = 2232.4 and 3080.0 
MeV are shown in Fig. 2. Good agreement between 
data and MC samples is observed, and a better reso¬ 
lution is achieved with increasing c.m. energy due to the 
smaller effects on the small angle multiple scattering. A 
c.m. energy dependent requirement, i.e., 9 p p > 178° at 
y/s < 2400.0 MeV, and 9 p p > 179° at /s > 2400.0 MeV, 
is further applied. Figure 3 shows the distribution of 
the momentum of proton or antiproton at c.m. energies 
y/s = 2232.4 and 3080.0 MeV. A momentum window of 
5 times the momentum resolution, |p mea — Pexpl < 5 a p , 
is applied to extract the signals, where p mea and Pexp are 
the measured and expected momentum of the proton or 


antiproton in the c.m. system, respectively, and a p is the 
corresponding resolution. 


B. Background study 

The potential background contamination can be classi¬ 
fied into two categories, the beam associated background 
and the physical background. 

The beam associated background includes interactions 
between the beam and the beam pipe, beam and resid¬ 
ual gas, and the Touschek effect [33]. Dedicated data 
samples with separated beams were collected with the 
BESIII detector at y/s = 2400.0 and 3400.0 MeV; these 
are used to study the beam associated background. Since 
the two beams do not interact with each other, all of the 
observed events are beam associated background, and 
can be used to evaluate the beam associated background 
at different c.m. energies by normalizing the data-taking 
time and efficiencies. No events from the separated beam 
data samples survive the signal selection criteria. Con¬ 
sidering that the normalization factor is less than 5 for 
most of energy points (other than 3.08 and 3.65 GeV), 
the beam associated background at all c.m. energy points 
is negligible. 

The physical background may come from the e + e - 
annihilation processes with two-body final states, e.g. 
Bhabha or di-muon events, where leptons are misiden- 
tified as protons or antiprotons, or processes with multi¬ 
body final states including pp, e.g. e + e~ -A pp 
The contamination from physical background is eval¬ 
uated by MC samples, and are listed in Table I for 
y/s = 2232.4 and 3080.0 MeV, respectively. 

The number of the surviving background events after 
normalization, is very small at the low c.m. en¬ 

ergies and can therefore be safely neglected. However, 
at higher c.m. energies {y/s > 3.40 GeV), due to the 
rapid decrease of the cross section of e + e _ -A pp, the 
background level which is mainly from Bhabha events is 
higher, and N//// needs to be corrected for. 

The ratio of pp invariant mass and the c.m. energy, 
Mpp/y/s, from data and MC has been compared and is 
shown in Fig. 4 at different c.m. energies. The integral 
luminosity of the data set at each c.m. energy is listed in 
Table II. There is good agreement between data and MC 
simulations. The signal yields are extracted by counting 
the number of events and are listed in Table II, where 
the quoted uncertainties are statistical only. The data 
sample at 3550.7 MeV is a combination of three data 
sub-samples with very close c.m. energies, y/s = 3542.4, 
3553.8, 3561.1 MeV, and the value of 3550.7 MeV is the 
average c.m. energy weighted with their luminosity val¬ 
ues. 
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FIG. 2. Opening angle distributions between proton and antiproton at the c.m. energies of (a) 2232.4 MeV, 
and (b) 3080.0 MeV. The dots with error bars are data, the histograms represent the distributions of signal MC 
samples. The arrows show the selection applied. 




p ( GeV/ C ) 

p/p 


FIG. 3. Momentum distribution of the proton or antiproton at the c.m. energies (a) 2232.4 MeV, and (b) 3080.0 
MeV, two entries per event. The dots with error bars are data, the histograms represent the distributions of signal 
MC samples. The arrows show the momentum window requirements. 


TABLE I. Physical background processes estimated from the MC samples at = 2232.4 and 3080.0 MeV. is the number 
of generated MC events, N™? is the number of events remaining after the selection criteria, a is the production cross section in 
the e + e _ annihilation process, which is obtained using the Babayaga generator for Bhabha, di-muon, and di-photon processes, 
and from the previous experimental results for others processes [34, 35]. and A^p are the estimated upper limit at 

the 90% confidence level (C.L.) and the normalized number of background events. 



a / S * 

2232.4 MeV (2.63 pb" 1 


a / s - 

3080.0 MeV (30.73 pb“ 

l ) 

Bkg. 


ct ( nb ) 

TyMC 

uplimit 

lyiViU 
- 1 ’nor 

Ngz ( xl 0 B ) N™ 

O' (" b ) ^unlimit 

lyiViU 

nor 

e^e - 

9.6 

0 

1435.01 

< 0.96 

0 

39.9 

1 

756.86 

< 2.54 

i 


0.7 

0 

17.41 

< 0.16 

0 

1.5 

0 

8.45 

< 0.42 

0 

77 

1.9 

0 

70.44 

< 0.24 

0 

4.5 

0 

37.05 

< 0.62 

0 

7 T~^7T~ 

0.1 

0 

0.17 

< 0.01 

0 

0.1 

0 

< 0.11 

< 0.02 

0 

K+K- 

0.1 

0 

0.14 

< 0.008 

0 

0.1 

0 

0.093 

< 0.02 

0 

pp7T° 

0.1 

0 

< 0.1 

< 0.006 

0 

0.1 

0 

< 0.1 

< 0.07 

0 

pp7r°7r° 

0.1 

0 

< 0.1 

< 0.006 

0 

0.1 

0 

< 0.1 

< 0.07 

0 

AA 

0.1 

0 

< 0.4 

< 0.02 

0 

0.1 

0 

0.002 

< 0.001 

0 
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FIG. 4. Comparison of M pp /yfs distributions at different c.m. energies for data (dots with error bars) and MC (histograms): 
(a) 2232.4, (b) 2400.0, (c) 2800.0, (d) 3050.0, (e) 3060.0, (f) 3080.0, (g) 3400.0, (h) 3500.0, (i) 3550.7, (j) 3600.2, (k) 3650.0, 
(1) 3671.0 MeV. The sample (i) is a combination of three data sub-samples with very close c.m. energies, y/s = 3542.4, 3553.8, 
3561.1 MeV, and the value of 3550.7 MeV is the average c.m. energy weighted with their luminosity values. 


C. Extraction of the Born cross section of 
e + e — —y pp and the effective FF 


The differential Born cross section of e + e pp can 
be written as a function of FFs, \Ge\ and \Gm\ [36], 


do"f iom ( S ) (X f~> C 

4s 


[|Gm(s)| 2 (1 + cos 2 0 P )+ 

4m 2 „ „ 

--|G_e(s)| 2 sin 2 d p ], 


( 4 ) 


where a « is the fine structure constant, /3 = 

/ 4m 2 

a/1 -— is the velocity of the proton in the e + e“ 

c.m. system, G = ^ i- C xp(—rra/p) ' s ^ le Coulomb cor¬ 
rection factor for a point-like proton, s is the square of 
the c.m. energy, and 9 P is the polar angle of the proton 
in the e + e _ c.m. system. We assume that the proton 
is point-like above the pp production threshold, meaning 
that the Coulomb force acts only on the already formed 


















































































































hadrons. At the energies we are considering here, the 
Coulomb correction factor can be safely assumed to be 
1. Furthermore, under the assumption of the effective 
FF |G| = \Ge\ = \Gm\ and by integrating over 9 P , it can 
be deduced: 


\G\ = 


O'Born 


86.83-1(1 + ^)’ 


( 5 ) 


where CTBorn is in nb and m p , s in GeV. 

Experimentally, the Born cross section of e + e _ — > pp 
is calculated by 


O'Born 


Nobs — Abkg 

L-e - (1 + 5)’ 


( 6 ) 


where N a b s is the observed number of candidate events, 
extracted by counting the number of signal events, iVbkg 
is the expected number of background events estimated 
by MC simulations, L is the integrated luminosity esti¬ 
mated with large-angle Bhabha events, e is the detec¬ 
tion efficiency determined from a MC sample generated 
using the Conexc generator [30], which includes radia¬ 
tive corrections (which will be discussed in detail in next 
paragraph), and (1 + 5) is the radiative correction factor 
which has also been determined using the Conexc gener¬ 
ator. The derived Born cross section o’Born, the effective 
FF |G|, as well as the related variables used to calculate 
o'Born are shown in Table II at different c.m. energies. In 
the table, the product value s' = ex (1+5) is presented to 
account for the effective efficiency. Comparisons of O’Born 
and |G| to the previous experimental measurements are 
shown in Fig. 5. Compared to the BaBar results [15], the 
precision of the Born cross section is improved by 30% for 
data sets with yfs < 3080.0 MeV, and the corresponding 
precision of effective FF is improved, too. 

From Eq. 4, it is obvious that the detection efficiency 
depends on the ratio of the electric and magnetic FFs, 
|Ge/Gm|, due to the different polar angle 9 p distribu¬ 
tions. In this analysis, the detection efficiency is evalu¬ 
ated with the MC samples. The ratio of \Ge/Gm\ is mea¬ 
sured for data samples at c.m. energies y/s = 2232.4 and 
2400.0 MeV, and for a combined data with sub-data sam¬ 
ples at = 3050.0, 3060.0, and 3080.0 MeV, which have 
close c.m. energy. The corresponding measured \Ge/Gm\ 
ratios are used as the inputs for MC generation. De¬ 
tails of the |Gb/Gm| ratio measurement can be found 
in Sec. HID. For other c.m. energy points, where the 
|Ge/Gm| ratios are not measured due to limited statis¬ 
tics, the detection efficiencies are obtained by averaging 
the efficiencies with setting |G_e| = 0 and \Gm\ = 0, 
respectively. The corresponding product values of detec¬ 
tion efficiencies and the radiative correction factors at 
different c.m. energies are listed in Table II. The interfer¬ 
ence of pp final states between e + e - annihilation and J/ip 
decay in the lower tail is assumed to be negligible [37]. 

Several sources of systematic uncertainties are consid¬ 
ered in the measurement of the Born cross sections and 


the corresponding effective FFs, including those of track¬ 
ing, PID, E/p requirement, background estimation, the¬ 
ory uncertainty from radiative corrections, FF model de¬ 
pendence and integrated luminosity. 


(a) Tracking and PID : The uncertainties of track¬ 
ing and PID efficiencies for proton/antiproton are in¬ 
vestigated using almost background-free control samples 
J/ip —> ppTr + Tr~ and ip(3686) — > tt + tt~ J/ip — > ■n + T:~pp. 
The differences of tracking and PID efficiencies between 
data and MC simulation is 1.0% per track, respectively, 
and they are taken as systematic uncertainties, (b) E/p 
requirement : The uncertainty of the E/p requirement 
is also estimated using the J/ip ppTT + n~ control sam¬ 
ple. The difference between data and MC in efficiency is 
found to be 1.0% applying the same E/p criteria on the 
proton sample, and is taken as a systematic uncertainty. 
(c) Background estimation : In the analysis, the back¬ 
ground contamination is estimated by the MC samples. 
An alternative method, 2-dimensional sidebands in the 
proton momentum versus antiproton momentum space, 
is applied to estimate the background contamination, and 
the difference is taken as the systematic uncertainty. The 
proton/antiproton momentum sideband region is defined 
by 6 a p < IPmea-Pexpl < U cr p , where p exp and a p are the 
expected momentum and resolution of proton/antiproton 
at a given c.m. energy, (d) Radiative correction : In 
the nominal results, the radiative correction factors are 
estimated with the Conexc generator. An alternative 
generator, Phokhara, is used to evaluate the theoret¬ 
ical calculation of the radiative correction factors, and 
the differences in the resulting products e' of detection 
efficiency and radiative correction factor are taken as the 
systematic uncertainty, (e) FFs model dependence : For 
those c.m. energies with measured \Ge/Gm\ ratios, the 
uncertainties on the detection efficiencies are estimated 
by varying the \Ge/Gm\ ratios with 1 standard deviation 
measured in this analysis. These systematic uncertain¬ 
ties are found to be less than 5.0%. For other c.m. en¬ 
ergy points, whose \Ge/Gm\ ratios are unknown, the 
uncertainties on the detection efficiencies are evaluated 
to be half of the differences between the detection effi¬ 
ciencies with setting \Ge\ = 0 or |Gm| = 0, respectively, 
which give larger uncertainties exceeding 10.0%. (f) In¬ 
tegrated luminosity : The integrated luminosity is mea¬ 
sured by analyzing large-angle Bhabha scattering pro¬ 
cess, and achieves 1.0% in precision. 


All systematic uncertainties are summarized in Ta¬ 
ble III. The total systematic uncertainty of the Born cross 
section is obtained by summing the individual contribu¬ 
tions in quadrature. The effective FF |G| is proportional 
to the square root of the Born cross section, and its sys¬ 
tematic uncertainty is half of that of the Born cross sec¬ 
tion. 
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TABLE II. Summary of the Born cross section <TBom, the effective FF |G|, and the related variables used to calculate the Born 
cross sections at the different c.m. energies \fs, where N 0 b s is the number of candidate events, IVbkg is the estimated background 
yield, s' = s x (1 ± 5) is the product of detection efficiency e and the radiative correction factor (1 + d), and L is the integrated 
luminosity. The first errors are statistical, and the second systematic. 


V3 (MeV) 

Aobs 

ATkg 

e' (%) L (pb- 1 ) 

O'Born (pb) 

|G| (xlO-^) 

2232.4 

614 ± 25 

1 

66.00 

2.63 

353.0 ± 14.3 ± 15.5 

16.10 ±0.32 ±0.35 

2400.0 

297 ± 17 

1 

65.79 

3.42 

132.7 ± 7.7 ±8.1 

10.07 ±0.29 ±0.31 

2800.0 

53 ± 7 

1 

65.08 

3.75 

21.3 ± 3.0 ± 2.8 

4.45 ±0.31 ±0.29 

3050.0 

91 ± 10 

2 

59.11 

14.90 

10.1 ± 1.1 ± 0.6 

3.29 ±0.17 ±0.09 

3060.0 

78 ±9 

2 

59.21 

15.06 

8.5 ± 1.0 ± 0.6 

3.03 ± 0.17 ± 0.10 

3080.0 

162 ± 13 

1 

58.97 

30.73 

8.9 ± 0.7 ± 0.5 

3.11 ±0.12 ±0.08 

3400.0 

2± 1 

0 

63.34 

1.73 

1.8 ± 1.3 ± 0.4 

1.54 ±0.55 ±0.18 

3500.0 

5± 2 

0 

63.70 

3.61 

2.2 ± 1.0 ± 0.6 

1.73 ±0.39 ±0.22 

3550.7 

24 ± 5 

1 

62.23 

18.15 

2.0 ± 0.4 ± 0.6 

1.67 ±0.17 ±0.23 

3600.2 

14 ±4 

1 

62.24 

9.55 

2.2 ± 0.6 ± 0.9 

1.78 ±0.25 ±0.35 

3650.0 

36 ±6 

4 

61.20 

48.82 

1.1 ± 0.2 ± 0.1 

1.26 ±0.11 ±0.07 

3671.0 

6± 2 

0 

51.17 

4.59 

2.2 ± 0.9 ± 0.8 

1.84 ± 0.37 ± 0.33 



M p p(GeV/c 2 ) 



M p p(GeV/c 2 ) 


FIG. 5. Comparison of (a) the Born cross section and (b) the effective FF |G| between this measurement and previous 
experiments, shown on a logarithmic scale for invariant pp masses from 2.20 to 3.70 GeV/c 2 . 


D. Extraction of the electromagnetic \Ge/ Gm\ 
ratio 


The distribution of the proton polar angle 9 p depends 
on the electric and magnetic FFs. The Eq. 4 can be 
rewritten as: 


F(cos9 p ) =iV norm [l + cos 2 0 P + 


4m| 

s 


f? 2 ( 1 


cos 2 Q p )\, 


( 7 ) 


where R = \Ge/Gm\ is the ratio of electric to magnetic 

FFs, and V norm = ?I^[IM+5M^R 2 }Gm{s) 2 is the 
overall normalization factor. Both R and N norm ( Gm{s )) 
can be extracted directly by fitting the cos 8 P distribu¬ 
tions with Eq. 7. The polar angular distributions cos 9 p 
are shown in Fig. 6 for y/s = 2232.4 and 2400.0 MeV, as 
well as for a combined data sample with sub-data sam¬ 
ples at yfs = 3050.0, 3060.0 and 3080.0 MeV. The dis¬ 
tributions are corrected with the detection efficiencies in 


different cos 9 p bins which are evaluated by MC simu¬ 
lation samples. The distributions are fitted with Eq. 7, 
and the fit results are also shown in Fig. 6. The fit re¬ 
sults as well as the corresponding qualities of fit, % 2 /ndf, 
are summarized in Table IV. The corresponding ratios 
R = \Ge/Gm\ are shown in Fig. 7, and the results from 
the previous experiments are also presented on the same 
plot for comparison. 

The systematic uncertainties of the \Ge/Gm \ ratio and 
\Gm\ measurements are mainly from background con¬ 
tamination, the difference of detection efficiency between 
data and MC, and the different fit range of cos 9 p . The 
small background contamination as listed in Table II is 
not considered in the nominal fit. An alternative fit with 
background subtraction is performed, where the back¬ 
ground contamination is estimated by the two-dimension 
sideband method, and the differences are considered as 
the systematic uncertainties related to background con¬ 
tamination. In the fit, the detection efficiency is eval¬ 
uated with the MC simulation. An alternative fit with 
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TABLE III. Summary of systematic uncertainties (in %) for the Born cross sections UBom and the effective form factor |G| 
measurements. 


77 ( MeV ) 

Trk . 

PID 

E/p 

Bkg . 

MC gen . 

Model 

Lum . 

Total (< r Born ) 

Total (| G |) 

2232.4 

2.0 

2.0 

1.0 

2.6 

0.4 

1.5 

1.0 

4.4 

2.2 

2400.0 

2.0 

2.0 

1.0 

2.0 

1.8 

4.5 

1.0 

6.1 

3.1 

2800.0 

2.0 

2.0 

1.0 

1.9 

7.5 

10.2 

1.0 

13.2 

6.6 

3050.0 

2.0 

2.0 

1.0 

2.2 

0.9 

4.0 

1.0 

5.6 

2.8 

3060.0 

2.0 

2.0 

1.0 

3.8 

0.1 

4.1 

1.0 

6.4 

3.2 

3080.0 

2.0 

2.0 

1.0 

0.0 

0.1 

4.3 

1.0 

5.3 

2.7 

3400.0 

2.0 

2.0 

1.0 

0.0 

7.8 

21.9 

1.0 

23.5 

11.8 

3500.0 

2.0 

2.0 

1.0 

20.0 

7.0 

12.9 

1.0 

25.0 

12.5 

3550.7 

2.0 

2.0 

1.0 

20.8 

9.0 

14.3 

1.0 

27.0 

13.5 

3600.2 

2.0 

2.0 

1.0 

35.7 

4.3 

11.6 

1.0 

37.9 

18.9 

3650.0 

2.0 

2.0 

1.0 

3.3 

0.9 

9.7 

1.0 

10.8 

5.4 

3671.0 

2.0 

2.0 

1.0 

33.3 

0.7 

13.3 

1.0 

36.0 

18.0 



COS0 p cos6 p cos6 p 


FIG. 6. Efficiency corrected distributions of cos 9 P and fit results for data at c.m. energies (a) 2232.4, (b) 2400.0 MeV and 
(c) a combined sample with c.m. energy at 3050.0, 3060.0 and 3080.0 MeV. The dots with error bars represent data. The solid 
line (black) represents the overall fit result. The dot-dashed line (in red) shows the contribution of the magnetic FF and the 
dashed line (in blue) of the electric FF. 


corrected detection efficiency which takes into account 
the differences in tracking, PID and E/p selection effi¬ 
ciency between data and MC is performed, and the re¬ 
sults in differences are taken as the systematic uncer¬ 
tainties. Fits with ranges [—0.8,0.6] and [—0.7, 0.7] in 
cos 6 p are performed, and the largest differences to the 
nominal values are taken as the uncertainties. Table V 
summarizes the related systematic uncertainties for the 
\Ge/Gm\ an d |Gm| measurements. The overall system¬ 
atic uncertainties are obtained by summing all the three 
systematic uncertainties in quadrature. 


(cos 2 6 P ) 


1— f 27TQ! pc cog 2 0^ + cog 2 0 p )\G M \ 
-‘’norm J 

4772 ^ 

—-(1 — cos 2 9 p )\R 2 \GM\ 2 ]dcos6 p . 


( 8 ) 


Calculating this within the interval [—0.8, 0.8] where 
the acceptance is non-zero and smooth, gives for the ac¬ 
ceptance correction: 


R = 


s (cos 2 6 P ) - 0.243 
4m 2 0.108- 0.648 (cos 2 9 V ) 


(9) 


and the corresponding uncertainty: 


As a crosscheck, a different method, named method of 
moments (MM) [38], is applied to extract the \Ge/Gm\ 
ratio, where the weighted factors in front of Ge and Gm 
may be used to evaluate the electric or magnetic FF from 
moments of the angular distribution directly. The ex¬ 
pectation value, or moment, of cos 2 9 p , for a distribution 
following Eq. 7 is given by: 


0.0741 s 

aR ~ R{ 0.167- (cos 2 9)Y 4m 2tT(cos2ep> ’ 
where <t < cos 2 8p ) is given by 


^(cos 2 6 P ) 



(cos 4 0p) — (cos 2 9 P ) 2 


( 10 ) 


( 11 ) 
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TABLE IV. Summary of the ratio of electric to magnetic FFs \Ge/Gm\, magnetic FF \Gm\ by fitting on the distribution of 
cos# p and method of moments at different c.m. energies. For the method of fitting on cos 9 P , the statistical and systematic 
uncertainties are quoted for \Ge/Gm\ and \Gm\, and the fitting quality \ 2 /n.d.o.f . is presented. Only statistical uncertainty 
is shown for the method of moments. 


Vi (MeV) 

\Ge/Gm\ 

\G m \ (xlO- 2 ) 

X 2 /ndf 



Fit on cos 6p 


2232.4 

0.87 ± 0.24 ± 0.05 

18.42 ± 5.09 ± 0.98 

1.04 

2400.0 

0.91 ±0.38 ±0.12 

11.30 ± 4.73 ± 1.53 

0.74 

(3050.0, 3080.0) 

0.95 ± 0.45 ± 0.21 

3.61 ± 1.71 ± 0.82 

0.61 


method of moments 


2232.4 

0.83 ± 0.24 

18.60 ±5.38 

- 

2400.0 

0.85 ± 0.37 

11.52 ± 5.01 

- 

(3050.0, 3080.0) 

0.88 ± 0.46 

3.34 ± 1.72 

- 



FIG. 7. The measured ratio of electric to magnetic FFs \Ge/Gm\ at different c.m. energy from BESIII (filled circles), BaBar 
at SLAC (open crosses) and PS170 at LEAR/CERN (open circles). 


In the analysis of experimental data, (cos 2 9 p ) and 
(cos 4 d p ) are the average of cos 2 9 p and cos 4 9 p which are 
calculated taking the detection efficiency event-by-event 
into account: 

_ i N 

(cos 2,4 9 P ) = cos 2 ' 4 9 P = — ^ cos 2,4 9 pi /£i, (12) 

V i=l 

where e, is the detection efficiency with the zth event’s 
kinematics as estimated by the MC simulation. 

The extracted \Ge/Gm\ ratios and \Gm\ by MM at 
different c.m. energies are also shown in Table IV, where 
\Gm\ is calculated by N norm in Eq. 7 using the mea¬ 
sured \Ge/Gm\ ratio. The results are well consistent 
with those extracted by fitting the distribution of polar 
angle cos 9 p , and the statistical uncertainty is found to 
be comparable between the two different methods due to 
the same number of events. 


IV. SUMMARY 


Using data at 12 c.m. energies between 2232.4 MeV 
and 3671.0 MeV collected with the BESIII detector, we 
measure the Born cross sections of e + e _ —> pp and ex¬ 
tract the corresponding effective FF |G| under the as¬ 
sumption \Ge\ = \Gm\- The results are in good agree¬ 
ment with previous experiments. The precision of the 
Born cross section with yfs < 3.08 GeV is between 
6.0% and 18.9% which is much improved comparing with 
the best precision of previous results (between 9.4% and 
26.9%) from BaBar experiment [15]; and the precision is 
comparable with those of previous results at yfs > 3.08 
GeV. The \Ge/Gm\ ratios and \Gm\ are extracted at the 
c.m. energies yfs = 2232.4 and 2400.0 MeV and a com¬ 
bined data sample with c.m. energy of 3050.0, 3060.0 and 
3080.0 MeV, with comparable uncertainties to previous 
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TABLE V. Summary of systematic uncertainties (in %) in the \Ge/Gm\ ratio and \Gm\ measurement. 


Source 

\Ge/Gm\ 

|g m | 

(MeV) 

2232.4 

2400.0 

(3050.0,3080.0) 

2232.4 

2400.0 (3050.0,3080.0) 

Background contamination 

1.1 

7.7 

3.2 

1.4 

7.7 

3.2 

Detection efficiency 

2.3 

1.1 

4.2 

2.3 

1.1 

4.2 

Fit range 

4.6 

11.0 

22.1 

4.6 

11.0 

22.1 

Total 

5.3 

13.5 

22.7 

5.3 

13.5 

22.7 


experiments. The measured \Ge/Gm\ ratios are close to 
unity which are consistent with those of the BaBar ex¬ 
periment in the same q 1 2 3 * 5 * 7 8 9 10 11 12 13 region. At present, the precision 
of the \Ge/Gm\ ratio is dominated by statistics. A MC 
simulation study shows that the precision can achieve 
10% or 3.0% if we have a factor of 5 or 50 times higher 
integrated luminosity. In the near future, a new scan 
at BEPCII with c.m. energy ranging between 2.0 GeV 
and 3.1 GeV is foreseen to improve the precision of the 
measurement on \Ge/Gm\ ratio in a wide range. 
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